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ABSTRACT 

We present A^-band (8-13 |4m) spectroscopic observations of the low-mass, embedded pre- 
main-sequence close binary system SVS13. Absorption features are clearly detected which 
are attributable to amorphous silicates, crystalline forsterite, crystalline enstatite and annealed 
Si02. Most intriguingly, a major component of the dust in the envelope or disc around SVS13 
appears to be SiC, required to model adequately both the total intensity and polarisation spec¬ 
tra. Silicon carbide is a species previously detected only in the spectra of C-rich evolved star 
atmospheres, wherein it is a dust condensate. It has not been unambiguously identified in the 
interstellar medium, and never before in a molecular cloud, let alone in close proximity to a 
forming star. Yet pre-Solar grains of SiC have been identified in meteorites, possibly suggest¬ 
ing an interesting parallel between SVS13 and our own Solar-System evolution. The unique¬ 
ness of the spectrum suggests that we are either catching SVS13 in a short-lived evolutionary 
phase and/or that there is something special about SVS13 itself that makes it rare amongst 
young stars. We speculate on the physical origin of the respective dust species and why they 
are all simultaneously present toward SVS13. Two scenarios are presented: a disc-instability- 
induced fragmentation, with subsequent localised heating and orbital evolution firstly anneal¬ 
ing initially amorphous silicates and then dispersing their crystalline products throughout a 
circumstellar disc; and a newly discovered shock-heating mechanism at the interface between 
the circumstellar and circumbinary discs providing the crystallisation process. One or both 
of these mechanisms acting on carbon-rich grain material can also feasibly produce the SiC 
signature. 

Key words: stars: formation - stars: individual (SVS13) - stars: variables: T Tauri, Herbig 
Ae/Be - ISM: dust. 


1 INTRODUCTION 


The Nebular Hypothesis, first formulated separately in the 18th 
century by Swedenborg, Kant and Laplace, has stood the test of 
time and is the generally favoured overall description of how the 
Solar System formed and evolved. In the young hot gaseous neb¬ 
ula, micron-sized dust grains are dragged by the gas, and when 
they collide they stick with one anot her and grow in size by form¬ 
ing fl uffy or porous aggregates (e. g.. IWeidenschillindUOOOl : IPopod 
2003|). Although the de tails are still being worked out (see e.g., 
Chiang & YoudiniboiOh . those aggregates continue to grow by col¬ 
lision until they become > 1 km in size (now called planetesimals), 
at which stage gravity takes over and fin ally ‘runaway’ growth 
of a few bodies leads to planet formation l lWeidenschilling||200d : 


* Based on data collected at Subaru Telescope, which is operated by the 
National Astronomical Observatory of Japan, 
t E-mail: tak@subaru.naoj.org 


IChiang & Youdinll2010h . It is therefore of great interest to study 
the mineralogy of circumstellar dust around young stars as it rep¬ 
resents the original constituents of planetesimals and hence of the 
rocky planets like our own Earth. 

SVS13 (star number 13 in IStrom jViba_& StromI [l976ll is 
a pre-main-sequence (PMS) star jAspinl 2003|) thought to be 
driving the Herbig-Haro objects 7-11 i Bachiller et al.l [l998l) 
near the reflection nebula NGC1333 in the Perseus molecular 
cloud. The bolometric luminosity is estimated to be ~52 Lq 
jMolinari. Liseau & Lorenzettilll993l : their estimate o f 115 Lp, at 
350 p c scaled to an updated distance of 235 ± 18 pc, iHorita et all 
120081) . Its optical and infrared brightnesses increased consider- 
ably sometime betw een 1988 December and 1990 September 
teisldffel et al.ll 19911) . The brightening was wavelength dependent, 
with the object becoming bluer as it underwent the change. Even 
after this period of relativ ely large increase in brightness (Amy ~ 
3.3 mag, Amk ~ 1.2 mag. lEisldffel et aLlll99lh . SVS13 showed a 
quasi-periodic near-infrared (NIR) brightness fluctuation with am- 
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plitude ~0.5 ma g and cycle ~500 days b etween 1990 October and 
1993 December lAspin & Sandelllll994ll . 

The 2.3-p.m CO overtone band-heads, indicative of the pres¬ 
ence o^Jiot_(~32000_^)_md_dense 10*® cm“®) molecular 

gas dCarr. Tokunaga & Naiitall2004h . probably in the shape of a 
circumstellar disc, have been observed in emission l Can! 1 19891 : 
lEislbffel et al.ll99ll:lBiscava et al.lll997l:ICarr et al.l200il . The CO 
band emission also displa yed variability on timescales as short as 
days teiscava et al.ll 199'^ . 

SVS13 is associated with a centimetre source (VLA 4, 
iRodriguez, Anglada & Curi^ll997h which has been resolved into 
a close binary (VL A 4A & VLA 4B) separated by ~0.3 arc 


_^_ b y 

sec, first at 3.6 cm (lAnglada. Rodriguez & TorrellesI [ 2 OOOI 1 then 
at 7 mm jAnglada et all 120041) . Although VLA 4A and 4B ex¬ 
hibit similar flux densities at centimetre wavelengths, VLA 4B is 
the dominant object (by a factor of more than 2) in the millime¬ 
tre range with its spectral index at these ra dio wavelengths mos t 
probably arising from thermal dust emission jAnglada et al.l2003) . 
The optical/NIR source SVS13 has been associated with VLA 4A 
by these authors based on older optical data; however, recently 
iLlodapp & Chinil j2014h identified SVS13 with VLA 4B, this time 
using the newer and more (both photometrically and astrometri- 
cally) accurate 2MASS catalogue (also, the epoch of the 2MASS 
observations happens to fall between the epochs of the two radio 
measurements that discovered and then confirmed the close bi¬ 
nary). The separation (~0.3 arcsec), at the assumed distance to the 
object of 235 pc, translates to 71 an, or slightly outside the Kuiper 
Belt in our Solar System. 

The first published mid-infrare d (MIR) photome t ry of 
SVS13 was conducted pre -outburst jCohen & Schwartj 1 19831 : 
iHarvev, Wilking & Joy| |l9 83). A decade later, without the kn owl- 
edge of the brightening, Liseau. Lorenzetti & Moling j 19921) ob¬ 
tained 10-[xm measurements but did not find the difference (from 
10 years before) substantial. The first A-band (8-13 pm) spec- 
tro scopy (in fact, spect ro-polarimetry) of the object was performed 
by lAitken et alj jl993h . using the UCL spectro-polarimeter on the 
3.8-m UKIRT with an etfective beam diameter of 4.3 arcsec and 
spectral resolution R ~ 40. They remarked that the polarisation 
spectrum of SVS13 is ‘unusual’, in that it peaked between 11 and 
12 pm, not at the characteristic wavelength of 10.2 pm for silicate 
grains. Nevertheless, they considered the polarisation to arise from 
dichro ic absorption, albeit by an unusual com position. Wright_et_^ 
j 19991) re-analysed the same data-set from lAitken et alj j 199311 
publishing the ‘u nusual’ spectrum for the first time (see also 
ISmith et aUbOOCh . They concluded that the double-trough-shaped 
intensity spectrum and the unique polarisation spectrum can only 
be satisfactorily modelled using inclusions of SiC in an amorphous 
silicate matrix (this point will be discussed in detail in 34.1.2.1L 

In this paper, we present the results of new A-band spectro¬ 
scopic observations of SVS13. The dust mineralogy in the circum¬ 
stellar environment is investigated by fitting various dust emissiv- 
ities to the spectrum ( 33.21 1. We also discuss dust components that 
comprise the best-fit model ( 34.11 1 and suggest their possible origins 

(pit. 


2 OBSERVATIONS AND DATA REDUCTION 

Table [T] summarises the observations of SVS13 and the standard 
star HD20644. All measure ments were made us ing the MIR imag¬ 
ing spectrometer COMICS jKataza et al.ll200d) at the Cassegrain 
focus of the 8.2-m Subaru Telescope on Mauna Kea, Hawaii. The 


A-band low spectral resolution (R ~ 250) mode of COMICS 
utilises two Raytheon 320x240 Si:As IBC arrays, one as a slit- 
viewer and the other spectrometer, and both are cooled by a Sum¬ 
itomo 4-K Gifford-McMahon-type cryo-cooler but usually operate 
at around 7~8 K because of the self-heating. We used a 2-pixel 
wide slit to achieve the spectral resolution of R ~ 250 in the A- 
band. The pixel scale of the spectrometer in the spatial dimension 
is 0.165 arcsec so the slit width translates to 0.33 arcsec, which is 
comparable to the diffraction-limited image size at 10 |xm at 8-m 
class telescopes. 

Data reduction was carried out using the ira^] data reduction 
and analysis package. The chop-subtracted 2-D spectra were flat- 
fielded and then geometry distortion corrected. Wavelength calibra¬ 
tion was achieved using numerous telluric emission lines present in 
the A-band. For flux calibration and telluric absorptio n correction, 
we a lso observed a Cohen standard star HD20644 jCohen et al.l 
I 1999 I) . We furthermore imaged both the target and standard in three 
different bands (To = 8.8, 11.7, and 12.4 pm) to rectify possible 
slit efficiency discrepancy. The images of SVS13 and HD20644, 
in terms of their point-spread functions, were very similar to each 
other at all three bands [i.e., SVS13 is unresolved. Indeed, when 
observing in the spectroscopy mode, the slit was set to run East- 
West (instrument PA = 90*’) in an unsuccessful attempt to resolve 
the close binary]. In fact, both sets of images of the target and the 
standard at longer wavelengths, i.e., at 11.7 and 12.4 pm, at which 
they are even less susceptible to seeing because of its T“5 wave¬ 
length dependence, show the first Airy disc, a clear sign of the 
diffraction-limited condition. We therefore conclude that only one 
star (VLA 4B) in the close binary system has any MIR emission 
associated with it. iHodapp & Chinil j20I4l) could also only find a 
single object in their laser-guided adaptive-optics assisted H- and 
R-band images. As such, in the final stages of the data reduction, 
three pixels in the spatial direction were simply summed up (i.e., 
3 X 0.165 ~ 0.5 arcsec) to improve the signal-to-noise ratio. 


3 RESULTS AND MODEL FITTING 
3.1 Spectrum 

Figure[T]shows the COMICS A-band spectrum of SVS13. Note that 
error-bars only represent the standard deviation in the sky back¬ 
ground and not the uncertainty in the absolute flux calibration, 
which would probably amount to about 10 per cent (most of which, 
approximately 4 to 7 per cent, originates from the uncertainty in the 
Cohen flux standard template at these wavelengths). 

The double-trough shape described in the Introduction is 
clearly apparent; however, at the spectral resolution R ~ 250, it re¬ 
markably becomes a triple-trough, with a hump at around 9.5 p.m. 
The first general dip is most likely caused by the silicate absorption, 
while the second, with its deepest p art at about 11 ixm, c ould be due 
to the presence of SiC, as argued bv lWright et al.l jl999l) . The high¬ 
est spectral resolution so far attained on the object by COMICS, 
compared to the relativ e ly low resolution of UCLS [R ~ 40, s ee 
fig.3(b) in I Wright et al.l j 19991) and fig.2.4 in I Smith et al.l j2000t) 1. 
displays finer details not previously seen, which will be utilised be¬ 
low for model fitting. 


* IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Table 1. COMICS observation log of SVS13 and the standard star HD20644. Integration is the total on-source integration time in seconds, and PA is position 
angle of slit in degrees, measured East of North. 


Date 

Object 

Band Integration 

PA 

Airmass 

(UT) 



(sec) 

(°) 




Spectroscopy 




2009/11/03 

SYS 13 

N 

994 

90 

1.224^2.125 


HD20644 

N 

65 

0 

1.072-^2.659 



Imaging 




2009/11/03 

SYS 13 

N8.8 

401 

- 

1.081^1.092 



N11.7 

401 

- 

1.102^1.113 



N12.4 

401 

- 

1.149^1.163 


HD20644 

N8.8 

22 

- 

1.059 



N11.7 

21 

- 

1.038 



N12.4 

21 

- 

1.047 



Figure 1. COMICS A-band spectrum of SVS13. Note that eiTor-bars only 
represent standard deviation in the sky background. 


3.2 Model fitting 

3.2.1 Dust species 

Selecting dust species can be subjective exercise but is constrained 
by those historically identified in astronomical spectra or as (pre- 
Solar) components of meteorites and interplanetary dust particles. 
We therefore restrict ourselves in using the following dust species 
and their combinations. 


(i) Amorphous (astr onomical) silicates (as represe nted by the 
Trapezium emissivity, iForrest. Gillett & SteinI Il975h . Although 
their specific mineralogy has been extremely difficult to identify, 
they are th e most abundant dust grains in th e interstellar medium 
(ISM; e.g.. iMolster. Waters & KempejbOl^ . The emissivity has 
been derived by d ividing the observe d A-band spectrum of the 


Trapezium reg ion l Fonestet_^ 1975tl by a 250-K Planck black- 


body function jOillettet all 1975 1. __ 

(ii) SiC (carbon star emissivity, lAitken et al.l ll97s il. It is a ma¬ 
jor co mponent in carbon-rich stellar ejecta (e.g., iTreffers & Coherj 
Il974h and is also found in meteorites as pre-Solar (pr obable extra- 
Solar-System) grains (e.g.. ISematowicz et alj 1987 ); however, it 
has n ever been observed in the ISM jWhittet. Dulev & Martini 
1 19901) . We note that, using a n Infrared Spac e Observatory (ISO) 
SWS spectrum, a later study l lMin et al.ll2007h has found some ev¬ 
idence of the existence of SiC grains in the A-band TSM’ absorp¬ 
tion feature in the line of sight towards the Galactic Centre (GC) 
source Sgr A*; however, the large beam size (14 x 20 arcsec") 


of ISO encompasses a number of asymptotic giant stars (i.e., the 
most probable source of meteoritic SiC grains, e.g., Zinnei| . 




identified in several catalogues j[BluimSellmcn_&^e2ovHl^9d; 
lOtt. Eckart & Genzelll 19991 ; IClenet et alj 200 li ; Blum et al. l2003h . 
We therefore suspect it may b e contaminated and s o may not repre¬ 
sent the pure ISM extinction. I Whittet et alj lll99d) stu died A-band 
spectra of 10 GC sources from Roche & AitkenI l l985h . which were 
obtained using the UCL spectrometer with a moderate beam diam¬ 
eter of 4.2 or 4.3 arcsec, thus those data should be less susceptible 
to contamination by source confusion. Regardless of its existence 
in the ISM, it has never been detected in the circumstellar environ¬ 
ment of young stars0 The emissivity has been derived by dividing 
the observed A-band spectru m of the carbon sta r Y-Tau by a 800-K 
Planck black-body function jAitken et aljl97lj) . 

(iii) Crysta l line silicates (aerosol forsterite and enstatite, 

iTamanai et al.l 1 20061) . The ISO ‘crystalline silicate revolution’ 
revealed they are frequently present in stars young and old 
djager et al.lll998t) . Now, Spitzer has both expanded and refined 
the crystalline silicate database (e.g.. lOlofsson et ^l2009l) ; how¬ 
ever, once again, they have never bee n identified in the diffuse 
ISM dKemper. Vriend & TielensI |2004|) [a caveat here is that the 
work of lMin_et_^ (2007 ) men tioned above used the same data¬ 
set from Kemper et alJ d2004h . i.e., the ISO SWS spectrum of 
the GC source Sgr A", which may not describe the true na¬ 
ture of the ISM dust grains]. The two most abundant crystalline 
silicate species are forsterite (Mg 2 Si 04 ) and enstatite (MgSiOs, 
iMolster et all l2010h . Commercially available crystalline silicate 
dust powders were dispersed into a nitrogen gas stream to create 
non-embedded, free-flying particle samples for infrared extinction 
measurements dTamanai et al.ll2006l) . _ 

(iv) Amorphous and annealed SiOi dpabian et al.l200ol) . It is one 
of the most abundant compounds in the Earth’s crust b ut its features 
have not been observed in the diffuse ISM (see e.g.. iLi & Draind 
l2002h . It has however been ide ntified in a number of protoplanetary 
discs around young stars (e.g.. lSargent et al]|2009ah . For the crys¬ 
talline component, we use the n ano-particles of pre cipitated silica 
annealed at 1,220 K for 5 hours l lFabian et alj200(]|) . 


^ An over-abundance of carbon has been reported in th e debris disc of 
the archetypal young star yS-Pictoris iRoberge et alj|2()()m . Although SiC 
is not explicitly referred to, it may be posited that one of the most likely 
compounds which could form in such a carbon-rich atmosphere is silicon 
carbide. 
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Figure 2. Normalised emissivities/absorption coefficients of the selected 
dust species (see text). 

We use either empirical (amorphous silicates and SiC) or 
experimental (crystalline silicates and amorphous/annealed SiOa) 
results on an ensemble of presumably differently shaped/sized 
grains. This approach is slightly varied from a number of other 
mineralogical studies which usually adopt some optical constants 
and shape distributions t o calculate absorption coefficients (e.g., 
iBohren & Huffma^l 198.4) . It is unrealistic to expect dust grains to 
have regular shapes and indeed a sample of interplanetary particles 
consi sts of irregularly shaped porous aggregates (see e.g. jBradlevI 
1200.4) . Also, SVS13 is polarised in the MIR (i.e., aligned aspherical 
grains must be present. IXitken et^ll993h . hence we consider our 
approach to be at least a more pragmatic one. Normalised emissiv¬ 
ities/absorption coefficients of the materials used in this study are 
shown in Figure[2] 

3.2.2 Model 

The fitting function can be expressed as 
Fa oc Ba{T) X 

where Ba(T) is the Planck black-body function at temperature T, 
a = T is ths total optical depth of k absorption components 
and T, is the optical depth of each component. In order to keep the 
fitting as simple as possible, we did not consider, for example, an 
extra black-body component, self-ab sorption of a dust species, etc. 

We used an iraf task specfit (E is to fit the emis¬ 

sivities/absorption coefficients. The fitting results are summarised 
in Table As the Planck black-body function is always present, 
we distinguish between each model with the number of absorp¬ 
tion components (k-comp, where k is the number of different dust 
species) and the names of specific absorptive ingredients added [am 
and cr stand for amorphous and crystalline, and Trap Trapezium 
(amorphous/astronomical silicates), fors/fo forsterite, enst/en en- 
statite, respectively. Every model is also given a unique alphabetical 
identification index (A-I), as listed at the beginning of each row]. 

The goodness of fit is given by the reduced chi-square, = 
IV, where v is the number of degrees of freedom. Since there is 
only one emissive component (a black-body) and the rest are al¬ 
ways absorptive, which all naturally come after the black-body, the 
order in which the absorptive components appear should not matter 
(i.e., they are commutative). We have indeed obtained exactly the 
same results (minimum and t values) even when the order in the 
component list was altered. It is interesting to note that when the 
Sic (plus Trapezium) emissivity is combined with only one of the 



Figure 3. Upper panel'. N-band spectrum of SVS13 overlaid with the best- 
fit model (solid line). Lower panel'. Fitting residual plot. 

crystalline silicate species (forsterite or enstatite), their best-fit pa¬ 
rameter values are rather uncertain [by more than 10 per cent; see 
models C (2-comp+cr fors) and D (2-comp-t-cr enst) in Table |3; 
however, when all three (SiC, crystalline forsterite and enstatite, 
plus Trapezium) are combined, their best-fit values become signif¬ 
icantly more certain [by < 5 per cent; see model E (2-comp-l-cr 
fo&en)]. This is most probably a manifestation of the fact that it 
takes all three emissivities/absorption coefficients to fit well the 
complex spectral shape longward of ~10 pm. 

The general trend is clearly such that the more components 
there are, the better the fit, as one might expect. However, it stops 
when both amorphous and annealed Si 02 are added to the mix. 
Were the trend to continue, one would have expected to obtain the 
smallest xi for ths combination that involves all the dust species 
considered in this study. Not only did xl not improve but also the 
contribution from amorphous Si 02 became insignificant (in addi¬ 
tion, its best-fit T is poorly constrained, with the associated uncer¬ 
tainty at nearly 30 per cent, while all the other values show rea¬ 
sonable certainty at well below 10 per cent). Combining these find¬ 
ings with the fact that adding annealed Si 02 alone gives a better 
fit than the amorphous counterpart in its place, we consider model 
G (4-comp-bcr Si02, i.e., a black-body, amorphous silicates, SiC, 
crystalline forsterite and enstatite, plus annealed Si 02 ) the best-fit 
model with the fewest number of components. Eigure[3 shows the 
A-band spectrum of SVS13 overlaid with the best-fit model (solid 
line). 

The model’s ability to reproduce fine structures, especially 
the longward of ~10 pm, is remarkable, considering the complex¬ 
ity of the spectrum. The only obvious shortcomings of the model 
prediction are: the slightly shallower slope towards the short-end 
of the A-band window; not reproducing well the shape of the 
8.3-pm absorption feature; the slightly deeper absorption at about 
9.1 pm; slightly shallower absorption at around 9.8 pm; not pre¬ 
dicting adequately the sharpness of the 10.1-pm absorption signa¬ 
ture, et cetera (see also the fitting residual plot in the lower panel 
of EigureO- However, it has been known for some time that, for 
laboratory measured absorption spectra of dust species, the method 
of preparation of dust samples is critically important. Eor exam- 
ple, hand-grinding and ball-milling will produce different results 
jKoike et alj|201ol) . Annealing processes (at different temperatures 
and for different durations) also affect the outcome and depending 
on the treatment, peak positions and full-widths at half-maximum 
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Table 2. specfit results of various model fitting (see text for details). In Column 1, alphabetical model IDs and names of various models discussed in the text 
are listed. The goodness of is given in Column 2. In Columns 3, the black-body temperature (in K) is given. In Columns 4—15, optical depths of each 

dust component and their associated uncertainties (given by specfit, converted to per cent) are listed. 



Model 

xl 

TiK) 

Trap 

err 

SiC 

eiT 

Crfo 

err 

Cr en 

en‘ 

Am Si02 

eiT 

Cr Si02 

err 

A 

l-comp 

30.0 

319 

1.11 

2 % 











B 

2-comp 

10.5 

298 

0.94 

1 % 

0.18 

4 % 









C 

2-comp+cr fors 

8.33 

300 

0.86 

1 % 

0.14 

13% 

0.08 

34 % 







D 

2-comp+cr enst 

6.51 

301 

0.66 

4 % 

0.18 

13% 



0.23 

13% 





E 

2-comp+cr fo&en 

3.74 

302 

0.55 

2 % 

0.15 

4% 

0.10 

5 % 

0.24 

1% 





F 

4-comp+am Si02 

1.95 

313 

0.56 

< 1 % 

0.17 

2% 

0.12 

4 % 

0.24 

2% 

0.15 

6% 



G 

4-comp+cr Si02 

1.67 

309 

0.60 

1 % 

0.18 

3% 

0.11 

5 % 

0.21 

5% 



0.16 

2 % 

H 

4-comp+am&cr Si02 

1.67 

310 

0.59 

< 1 % 

0.18 

3% 

0.11 

6% 

0.21 

3% 

0.02 

29% 

0.14 

4 % 

I 

G-SiC 

14.7 

314 

0.56 

< 1 % 



0.18 

4% 

0.23 

4% 






of du st featur es will differ, some times significantly jKoike et al.l 
I 2 OI 0 I : see also iFabian et alJboOOh . When these further complica¬ 
tions are taken into account, it is even more impressive how good a 
fit the selected model components provided, considering the com¬ 
plexity of the features modelled and the limited number of dust 
species required. In the following discussion, we assume that the 
dust components present in the best-fit model exist in the circum- 
stellar environment of the PMS star SVS13. 


4 DISCUSSION 

4.1 Dust species (revisited) 

4.1.1 Amorphous (astronomical) silicates 

As mentione d earlier, they are the most abundant dust species in the 
diffuse ISM l lMolster et al.l 1201(1) . The broad, smooth absorption 
feature arising from the Si-O stretching mode, peaking at about 
9.7 pm, is the str ongest and be st studied infrared feature in the 
diffuse ISM (e.g.. IWhittetll200^ . although identifying their spe¬ 
cific mineralogy or composition has traditionally been extremely 
difficult, especially if accompanying 20-pm data is not available. 
Most, if not all, astronomical mineralogy studies contain amor- 
phous silicates in o ne form or another (see, for example, tab.l in 
iMolster et al.ll201(lh . Therefore, it is not surprising that the amor¬ 
phous silicate absorption appears prominently in the N-band spec¬ 
trum of a young star. However, we do not attempt to associate it 
with a particular mineral, but rather use a generic representation 
of amorphous silicates in a star-forming region as provided by the 
Trapezium emissivity. 


4.1.2 SiC 


Silicon carbide has never been observed in the diffuse ISM 
dWhittet et al.l Il990l) . although it is produced in carbon star 
ejecta and is undoubtedly injected into it. So where does it 
go? One simple explanation may be that dust grains are de- 
stroyed in the ISM by, for exampl e, supernova shock waves (e.g., 
Ilones, Tielens & Hollenbai^ 1996h. If the sub stantial gap between 
the dust lifetime (~1Q^ vr. llones et aDll994 1996) and the dust 
injection timescale (from dust formation in the outflows of dying 
stars to being incorpo r ated into new star-form ing regions; ~10® yr, 
iDwek & Scald 1 19801 ; Ilones & Tielend 1 19941) is valid, then SiC 
grains may also have been destroyed, hence no 11.3-pm SiC fea¬ 
ture detected in the ISM. However, as mentioned earlier, pre-Solar 


SiC grains (and many othe rs) have been found (e.g. , in a carbona¬ 
ceous chondrite Murray bv iBematowicz et alJll987l) so some dust 
grains must survive the harsh ISM. 

SiC has in fact also been suggested as a dust component in 
a number of comets dOrofino, Blanco & Fonti|[l994l) . Good fits to 
the A-band spectra of three comets (IP/Halley, C/1987 PI Brad- 
field, and 19861 Wilson), which all exhibited two broad peaks cen¬ 
tred at about 9.7 and 11.3 qm| were obtained using a laboratory 
measured spectrum of synthetic amorphous olivine (i.e., silicate, 
Blanco et aljl99ll) and the complex refractive index of a-SiC from 
Pegourid ( Il988l) . It should be noted, however, that the 11.3-p.m 
feature in comets has generall y been attributed to cry stalline sili¬ 
cates (e.g., for Halley’s comet. ICampins & Rvan|[l989l) . While the 
available data is limited to the A-band spectrum alone, its iden¬ 
tification unfortunately remains rather ambiguous as their signa¬ 
tures are often overwhelmed by that of amorphous silicates. Only 
when measurements at other wavelengths (especially at 20 qm and 
longer) are on hand, the crystalline silicate assignment becomes 
more robust (e . g., fo r Comet Hale-Bopp using the ISO SWS data, 
ICrovisier et al.lll99'^ ; 1997). Those three comets (and Hale-Bopp 
also) are all long-period Oort Cloud comets which are thought to 
have originated just beyond the ice giant formatio n region (> 15 au, 
IChamoz & Morbidellil200^ ; lLevison et al.ll2008h . 

The most conclusive evidence for the presence of pre-Solar 
Sic grains in comets came from dust samples, captured and re- 
tumed to Earth by NASA’s Stardust mission to Comet 81P/Wild 2 
dMessenger et al.l l2009l) . which belongs to short-period lupiter 
Family of comets (IFCs), which are likely to have formed in the 
Kuiper Belt (> 30 au). The discovered 300-nm SiC grain is unique, 
in that it is the only pre-Solar grain so far found in the silica aerogel 
dust collector, possibly reflecti ng the complex p reparation required 
for the aerogel track analysis dFloss et alj|2013l) . and no pre-Solar 
SiC grains have yet been detected in the impact craters on the sur¬ 
face of exposed aluminium foils. The abundance of pr e-Solar SiC 
grain s in 81 P/Wild 2 has been estimated to be ~45 ppm dFloss et al.l 
l2013h . which is consistent with that in insoluble organic matter 
residues fro m a number of different classes of chondritic meteorites 
(10-55 ppm. lDavidson et al]l2009l) . These findings reveal that SiC 
grains must have been available in a wide range of distances from 
the Sun, from the carbonaceous chondrite formation region of the 
inner few au to the comet-forming outer reaches of the proto-Solar 
nebula. 


^ Comet Wilson showed an extra, unidentified e mission feature at aro und 
12.5 qm but it was omitted from fitting procedure dOrofino et al.lll99^ . 
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Interestingly, pre-Solar SiC grains found in another carbona¬ 
ceous chondritic met eorite Murchison, for ex ample, showed little 
evidence of scarring jBematowicz et al]|2003h due, for instance, to 
grain-grain collisions. Were the inters tellar environmen t as destruc¬ 
tive as it has been described (e.g., bv llones et al.lll99^ . such hall¬ 
marks should certainly be present. Then again, whether a dust grain 
survives is probably largely statistical in nature and it is perhaps 
only natural that those th at do show little damage. Nonetheless, 
when llones & NuthldTOllh re-examined the dust lifetime using re¬ 
vised uncertainty estimates, they found it to be comparable with the 
dust injection timescale. 

IWhittet et alj ( Il990h proposed an alternative solution to the 
missing SiC problem: oxidation. They argued that the ejected SiC 
grains could be selectively destroyed by surface oxidation in the 
0-rich interstellar environment. Indeed, partial oxide layers on 
some pristine pre-S olar SiC grains from the Murchison meteo rite 
have been reported dCroat, Stadermann & Bematowical2010ah . In 
volatilisation experim ents of SiC grains unde r Solar nebula-like 
oxidising conditions jMendvbaev et al.l l2002h . either continuous 
or partial Si 02 layers are formed, depending on oxygen fugaci- 
ties (continuous at higher and partial at lower fugacities, or more 
and less oxidising, respectively). Similarly, though probably un¬ 
intentionally, an oxide layer had grown on the surface of SiC 
nano-particles after they were taken out of a vacuum chamber 
Jciement et al.l |200^ . However, even though a Si02 layer does 
add new features due to the presence of the mantle material, it 
does not appear to suppress very much the 11.3-ttm SiC fe ature 
JPosch, Mutschke & Andersenll2004l : IZhang. Jiang & Lill2009h . 

In a similar study to ICroat et ^ ( 1201 Oah . 

l 2009ah and 

identified coatings 


Croat . Bernatowicz & Stad ermannI _ 

Croat. Lebsack & Bematowiczl ( 1201 Och 


mainly consisting of carbonaceous material on the surface of some 
pristine pre-Solar SiC grains also from Murchison. Furthermore, 
once again from the same meteorite, some graphite spherules 
were found to contain SiC grains ( IBernatowicz et al.l 1 1996 


Croal & SladermannI 20061; Hvnes 

Croat & Bernatowicz 

2007; 

Croal. Sladermann & Bernatowicz 

20081; ICroat & Stadermann 


SiC-containing graphites have also bee n found i n sam¬ 
ples from another meteorite, O rgueil dCroat et alj l2009bl : 
ICroat. Bernatowicz & Jadhi^ |2014|) . The discoveries of these 
SiC-containing graphites are quite significant as they are clear 
evidence that pre-Solar SiC grains either formed earlier than or at 
about the same time as the graphites. One spherule in particular 
(only one of a kind so far), discovered by ICroat & Staderma^ 
(120081) . houses a SiC grain at its centre, most probably indicating 
the central SiC acted as a nucleation core (i.e., SiC may have been 
formed first). It should be noted, however, that the number of SiC 
inclusions hitherto found is rather limited; only about a dozen 
graphite spherules that contain SiC grains have been reported (see 
references listed ab ove), whilst grap hites with internal carbides 
(predominantly TiC. ICroat et alj2008l) number ~80 to date, despite 
the fact that the abundance of Ti is very much lower than that of Si. 

These results (the scarcity of SiC-containing graphites and 
much more numerous TiC inclusions) can be explained by a rel¬ 
ative condensation sequence of solids in carbon star outflows; 
in conditions typically found for carbon stars, the condensa- 
tion sequence is though t to be TiC, graphite, and then SiC 
dLodders & FeglevI Il995h . There have however been some sug¬ 
gestions that, under certain circums t ances, SiC is formed be fore 
graphite dChigai & Yamamot3l2003l ; Ivasuda & Kozasair2012h . In 
fact, condensation in this order (SiC before graphite) has been in¬ 


voked to acc ount for the lack of the 11.3-trm SiC feature in the 
diffus e ISM dFrenklach. Carmer & Feigelsord 1 19891 ; iKozasa et al.l 
1 19961) . As SiC grains formed in the vicinity of carbon-rich stars 
travel further away and become cooler, a carbon layer develops 
on the grains and hides the SiC feature. Su ch scenario was in 
fact briefly mentioned by IWhittet et alj dl99C)h . and unlike the ox¬ 
ide coating discussed earlier, the car bon mantle seems to suppress 
succe ssfully the 11.3-p.m signature dKozasa et alJll993 ; IPapoul^ 

l2008h . 

A short letter has been published in a surface science 
journal that reports the dissolution of the amorphous carbon 
layer into the SiC core, when the core-mantle grain was heated 
to 600°C (873.15 K) und er high vacuum (3 x 10“* Torr, 
iKimura. Saito & Kaitd 200^ . At 800°C (1,073.15 K), the surface 
layer completely disappeared without altering the lattice structure 
of the particle. Somewhat surprisingly, they have further observed 
that when the SiC grain was cooled back to room temperature, 
still under high vacuum and thus without any significant contam¬ 
inants, the a morphous carbon coatin g re-emerged (i.e., the process 
is reversible, Iy. Kimura et alj|2003l) . The temperature reached in 
their experiment (1,073.15 K) is above the glass transition temper¬ 
ature of a morphous analogues of enstatite (MgSOs) composition 
(1,040 K. iRoskosz et al.ll201 ih . As crystalline dust species (crys¬ 
talline forsterite and enstatite, and annealed Si02; see 34.1.31 & 
14.1.4l l are also present, temperatures reached in the circumstellar 
environment of SVS13 must have been high enough to produce 
those and remove the carbon coating from the SiC grains. Under 
presumably much more severe conditions encountered, the amor¬ 
phous carbon layer could still be absent from the Si C core even 
at suffi ciently low temperatures (< 100 K at > 15 au. IPoteet et al.l 
l201lh for the grains to exhibit the SiC absorption feature. Recall 


Wright et alj d 19991) concluded that SiC grains are most likely 


that 

to be included in an amorphous silicate matrix (see also 34.1.2.1k 
To our knowledge, it has never been established whether the 
re-heated SiC particles would exhibit (in absorption) exactly the 
same feature as that seen towards carbon stars. In any event, it 
is known to vary somewhat from one carbon star to another; so 
much so that. lSneck. Thompson & Hofmeisted ( l2005h have coined 
a term ‘~ll-p.m feature’ (i.e., with the approximation mark 
and fewer significant figures) to emphasise its variability. Also, a 
self-absor bed SiC band has been detected towards extreme car¬ 
bon stars ( Speck, Barlow & Skinnej|l997l ; ljusttanont et al.|[l997l ; 


iHonv. Waters & Tielensll2002h . where amorphous carbon might be 

expected to dominate. Therefore, although this scenario (resur¬ 
gence of the Sic feature in a warm environment) may plausibly 
explain the situation in SVS13, it should be treated with caution. 

However, we highlight the fact here that if we take out the 
SiC component from the best-fit model (I = G - SiC = a black- 
body, amorphous silicates, crystalline forsterite and enstatite)3 the 
resulting fitting value (14.7) is better than that for the most basic 
model A (1-comp = a black-body plus amorphous silicates; 30.0) 
but worse even than the one for the simple model B (2-comp = 
a black-body, amorphous silicates, and SiC; 10.5), signifying need 
for the SiC inclusion. See Figure|4]which compares models B and I. 
It is clear, apparent in the residual plot in the lower panel, that the fit 


'* When SiC was taken out of the best-fit model, the contribution from an¬ 
nealed Si 02 became rather insignificant or even negative (i.e., turned emis¬ 
sive). We have therefore decided to exclude it from the fitting procedure 
in this test case. The inference drawn in this section would still be valid 
regardless of the crystalline Si 02 inclusion. 
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Figure 4. Upper panel'. W-band spectrum of SVS13 overlaid with models 
B and I. Lower panel: Fitting residual plot for the two models. 


for model I becomes much worse than that for model B longward of 
~10 ttm, especially around 11.3 p.m, where the contribution from 
the SiC emissivity is significant. 

4.1.2.1 Polarimetry perspective Similar to the conventional 
sp ectrum, the MIR po laris ation spectrum of S VS13, first discussed 
in lWright et al.l lll999h and lSmith et alJ (l2000h . is perhaps the most 
unique ever observed (see Figures |5] and [^. Given the peak occurs 
beyond 11 pm and the overall profile has an approximate ‘tilde’ 
shape, proper t ies similar to those predicted for polarised emission 
l lMartirj|l975l: lAitkenlll989tl . at first glance it is not even obvious 
whether the spectrum is representative of dichroic absorption or in¬ 
stead polarised emission. Also, that the total flux density spectrum 
is obviously dominated by absorption does not necessarily imply 
the polarisation profile would also be indicative of absorption. In¬ 
stances can occur where the polarisation is dominated by emission 
but the conventio nal flux density s hows deep absorption, an effect 
first predicted by iThronso an d seen in a few so urces (e.g., 

NGC7538 IRSl and RCW57 IRSl in ISmith et alboooh . 

Fortunately in the case of SVS13 there is additional polarisa¬ 
tion data in the NIR, up to almost 3.5 pm, which can be compared 
with the MIR jHollowav et al.ll2002l and references therein). Most 
specifically the measured polarisation position angles in the two 
spectral regions are almost identical at ~50 ± 5 degrees. The NIR 
polarisation is almost certainly due to dichroic absorption since it 
rises through the 3-pm water-ice band. There can be no polarised 
emission through this band because the ice would obviously melt at 
the 500-1,000-K temperatures characteristic of NIR emission (and 
it is unclear if relevant grain alignment mechanisms could operate 
at such a high temperature). Thus, given the equality of the posi¬ 
tion angles, it is very likely that the MIR polarisation is also due to 
dichroic absorption. 

Having established the polarisation process as dichroic ab¬ 
sorption it remains to determine what type of dust can produce such 
a unique 8-13-pm profile, and especially a peak beyond 11 pm. 
Possible candidate carriers with bands around 11 pm and longer, 
inferred in other astronomical spectra and/or found within pre- 
Solar meteoritic grains, include crystalline silicates and SiC, as 
shown in Figure [T] as well as water ice, Polycyclic Aromatic Hy¬ 
drocarbons (PAHs), aluminium oxides and carbonates. Polarisation 
models have been run for all of these potential constituents, either 


as inclusions mixed in with an amorphous silicate host matrix or as 
a confocal mantle coating an amorphous silicate core. 

Details of the modelling will be presented elsewhere (Wright 
et al. in preparation), but are based on the Rayleigh approxi- 
m ation for spheroidal particles . At M IR wavelengths, according 
to ISomsikov & Voshchinniko^ ( Il999h . this is adequate for grain 
radii (of an equivalent volume sphere) up to at least 0.5 pm, 
where the relative difference between exact and Rayleigh solu¬ 
tions is less than 10 per cent for the water ice coated grains of 
their work. Relevant formulae for the absorption and polarisa- 
tio n cross sections of b oth b are and mantle d spher oids are given 
in IPraine & Led j 19841) and iLee & Draind ( Il985ll . whilst those 
for an ‘average’ dielectric function of a mixture, using effective 
medium theory such as the Maxwell-Gamett rule, are provided in 
iBohren & Huffma^ ( Il983h . In all cases we have used an oblate 
core with principal axis ratio of 2, paramet ers which are reasonably 
constrained from other observation s (e.e.. lHildebrand & DragovanI 
ll995l:lDraine & Allaf-Akbarill2006h . 


The amorphous silicate optical constants are those of 
‘astron omical silicate’, first formula ted by Draine_&J^ej 
(Il984l). and subsequ e ntly updat ed by iLaor & Drainel (Il993h . 
IWeingartner & Draiiiel (1200ll) and IPraind ( l2003h . The MIR com¬ 
ponent of these optical data sets was constructed to reproduce 
the spectrum of Trapezium, so our choice is consistent with the 
emissivity function we used to model the conventional (total flux 
density) spectrum in 33.21 

As might be expected based on the absorption coefficients 
shown in Figure[2 polarisation features of crystalline silicates, suc h 
as inclusions of the crystalline olivine of iMukai & Koikd ( Il99(lh . 
are simply too narrow to match the observed SVS13 polarisation 
spectrum. They instead produce a spectrum which looks much like 
the polarisation profile of t he massive embedded young stellar ob- 
jects (YSOs) AFGL259 1 jAitken et al.lll988l: IWright et al.lll999l) 
and IRAS 13481-6124 jWright et alj l2008h . two sources with a 
sharp polarisation feature at ~11.2 pm superposed on the much 
broader amorphous silicate band. Similar comments hold for PAHs 
and carbonates. Whilst its resonance is relatively broad, aluminium 
oxide inclusions merely produce a shoulder, or lon g-wavelength 
wing , on the amorphous silicate polarisation profile dWright et al.l 

l2ooii . 


The long-sought but rarely identified 12-pm librational band 
of water ice, occurring between ~12 and 13 pm for the c rystalline 
and amorphous phases respectively dMaldoni et al.|[l993) . may be 
thought to be a strong candidate. SVS13 does have the character- 
istic 3.1-pm H 2 O ice absorption band, which Ip^ise et ah 1 120031) 
note is more like that from crystalline than amorphous ice. Its op¬ 
tical depth however is not particularly high, only around 0.55. This 
compare s to values of around 1 to 3 . 5 in a sample of embedded 
YSOs in ISmith. Sellgren & Tokuna^ dl989l) . Their study included 
another source with an unusual 8-13-pm spectrum, namely the 
eastern component of the double source AFGL961, with a 3.1-pm 
optical depth of 2.46. 

Indeed, the only previous identi fication of th e librational band 
was made for AFGL961, firstly by ICoxI d 19891) using IRAS and 
then bv ISmith & Wright d201lh for e ach component using Spitzer. 
This h as recently been questioned bv iRobinson. Smith & Maldonil 
d2012h on the basis of radiative transfer models suggesting the 
8-13-pm silicate feature could be self-absorbed, leading to a lo¬ 
cal minimum between the 10- and 20-pm silicate bands occurring 
around 13 pm and thus mimicking ice absorption. However, their 
models only do an average job of reproducing the overall shape of 
the observed 2-20-pm spectrum and profile of the 8-13-pm fea- 
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Figure 5. SVS13 polaiisation spectrum (solid circles) overlaid with vaiious polarisation models (solid lines) calculated for spheroids in the Rayleigh ap¬ 
proximation. The data around 9.5 ]im are compromised by telluric ozone. Insets show the optical depth profile extracted from the conventional flux density 
spectrum. The dashed line in all panels (insets) corresponds to the polarisation (absorption) expected from ‘astronomical silicate’ of|D rain 3 <2003h . In (b) and 
(d) are models with a mantle of amorphous (b) and crystalline (d) water ice, and a mantle-to-core volume ratio of 0.5. Models in (a) and (c) are for silicate-ice 
mixtures, again using amorphous (a) and crystalline (c) ice with a volume fraction of 0.375 and a Maxwell-Garnett mixing rule for the dielectric function of 
the effective medium. In all cases the grain cores are oblate with an axial ratio of 2:1. 


ture. Also, they only modelled AFGL961E and not AFGL961W, 
which has a prominent silicate emission feature but still shows the 
same apparent absorption at 13 p.m (as well as CO 2 ice absorptio n 
at 15 pm). Whatever the case, as presented in iBoogert et al.l ( l2008h . 
there is sufficient evidence that the librational band has been de¬ 
tected toward at least a few YSOs (although none look much like 
SVS13 itself). 

Thus, polarisation models with water ice either as a mantle 
or as inclusions were attempted. Ice existing as a mantle on sil¬ 
icate cores is the ‘sta ndard’ model for dust in molecular clouds 
(e.g.. lGibb et al.ll2004l) . On the other hand, the idea behind ice as 
inclusions is that if silicate grains are initially porous or fluffy, then 
within a cold molecular cloud (or disc) water ice may condense 
to fill the pores. Alternatively, ice-mantled grains may collide and 
stick together, resulting in a similar type of grain structure. 

Whatever the case, core-mantle grains cannot produce the po¬ 
sition of the polarisation peak in SVS13. See Figure|5] which shows 
in (b) and (d) models using a mantle of respe ctively amorphous and 
crystalline water ice (optical co nstants from iLeger et all 19831 and 
iBertie, Labbe & Whalleviri969L respectively). The mantle-to-core 


volume ratio is 0.5, but in fact mantles do not sufficiently change 
the overall polarisation profile even up to mantle-to-core volume 
ratios of 4 or more. The 10.3-pm peak of bare amorphous silicates 
is still clear, and indeed sharper and higher than the water ice peak 
around 12 pm, whilst the observed SVS 13 peak lies in the trough 
between them. 

On the other hand, as seen in Figure (a) and (c), a 
silicate-ice mixture can reasonably match at least the overall 
shape and peak position of the SVS 13 polarisation spectrum for 
an ice-to-silicate volume ratio of 0.375. A serious flaw how¬ 
ever is that it predicts large polarisation in the 3-pm water 
ice band, e.g., an excess of ~6-8 per cent over the contin- 
uu m compared t o an observed value of < 1 .5 per cent described 
in lAitkenI l ll996h and IChrvsostomou et ^ ( Il996l) . This assumes 
that grains are small enough that the Rayleigh approximation re¬ 
mains valid at 3 pm, reasonable for grai n sizes < 0.5 pm. See 
tab.2 in ISomsikov & Voshchinnik^ ( Il999h . where the relative dif¬ 
ference between exact and Rayleigh solutions only increases to 
10-15 per cent at 3 pm for either oblate or prolate grains with an 
‘equivalent volume’ spherical radius of 0.5 pm. 
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Figure 6. SVS13 polarisation spectrum (solid circles) overlaid with various polarisation models (solid lines) calculated for spheroids in the Rayleigh approxi¬ 
mation. The data ai'ound 9.5 p-m are compromised by telluric ozone. Insets show the optical depth pro file extracted f rom the conventional flux density spectrum. 
The dashed line in all panels (insets) corre sponds to the polarisation (absorption) expected from the lPrain^ j2003h ‘astronomical silicate’. Panels (a)-(c) use 
SiC optical data from Pitman et alJ ^2008ll for cubic (a), h exagonal e-ray (b) and hexagonal o-ray (c) sam ples. Panel (d) u ses a-SiC refractive indices from 
IChovke & PalikI jl985 ), panel (e) ‘astronoiuical SiC’ from iLaor & Drain3119931) and panel (f) cr-SiC from Pegourid il988l) . In all cases the grains are oblate 


with an axial ratio of 2:1. 


Included in each panel of Figurej^is an inset showing the opti- 
ca l depth extracted fr om the UCLS spectrum of S VS 13 as presented 
in ISmith et al.l ( l200(lh . compared to the absorption cross section of 
the respective models. The optical depth was calculated by fitting a 
Planck function between the short and long wavelength ends of the 
spectrum, after correcting them for a grey component of emissiv- 
ity taken to be the Trapezium emissivities at 7.5 and 13.5 pm. The 
inferred colour temperature is around 350 K, not too dissimilar to 
the model-based temperatures listed in Table|2] In no case does wa¬ 
ter ice, either as a mantle or as inclusions, match the optical depth 
profile. We therefore conclude that water ice is highly unlikely to 
be responsible for the unique SVS13 MIR spectrum, unless its op¬ 
tical properties are much different to those used here (which are 
extremely similar to those presented in many other publications). 

Thus, the only remaining alternative is SiC and indeed this 
provides the best match to the observed SVS13 MIR polarisa¬ 
tion spectrum. See Figure which shows polarisation models of 
SiC inclusions within an amorphous silicate matrix, using the 
Maxwell-Garnett mixing rule to calculate the effective dielectric 
function. Several different dielectric functions of SiC were trialed. 


namely: cubic (3C or j8) and two samples of hexagonal (6H or 
a) for E both perpe ndicular and parallel to the c-axis, derived by 
IPitman e t alj )l20 08h from single -crystal reflectance spectroscopy; 
Qf-SiC in Chovke & PalikI ( Il985h . also from reflec t ivity m easure- 
ments; an ‘astronomical SiC’ from iLaor & Draiii^ l ll993h derived 
from several different laboratory data sets; or-SiC determi ned from 
particle transmission measurements and presented by IPegouriel 


All can provide good matches to the overall profile of Ihe ob¬ 
served spectrum, with respective volume fractions of the SiC in¬ 
clusions of 0.1 for the Pitman et al. and Choyke & Palik data sets, 
and 0.25 for the Laor & Draine and Pegourie optical data. Whilst 
the models are all qualitatively similar, apart perhaps from the Pe¬ 
gourie case, the best match appears to be for the 6H (hexagonal) 
e-ray optical data of Pitman et al. (2008; Figure[^b). Specifically, 
the wavelength of peak polarisation and the profile shapes of both 
the polarisation and optical depth more closely follow the observa¬ 
tions. 

That the SiC polarisation signature dominates over that of 
silicates is testament to its stronger refractive indices and band 
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10 Fujiyoshi et al. 


strength, i.e., absorption cross section per unit volume. In at least 
the crystalline samples this is ac centuated by a polarisation rever¬ 
sal ac ross the strong resonance iMartinll 19751 : iHong & Greenberd 
Il978h . which essentially subtracts polarisation from the amorphous 
silicate at 10 p.m. 

In the interests of completeness, we computed models with 
Sic as a confocal mantle on a silicate core. These did not pro¬ 
vide anywhere near as good a match to the observed optical 
depth and polarisation profiles. Similarly, model s using the op¬ 
tical d ata of amorphous thin film s of SiC from iMutschke et al.l 
(Il999h and lLarruQuert et ^ 1201 ih did not match the spectra, be¬ 
ing either too broad and/or peaking at the wrong wavelength. 
Given their weaker refractive indices than for the crystalline phases 
they also necessitated a larger relative abundance of SiC (com¬ 
pared to silicates). Finally, we replaced the ‘astronomical silicate’ 
with a magnesium-iro n laboratory olivine, MgegFe, 2 Si 04 , from 
IPorschner et al.l ll995h . motivated by its good match to the di ffuse 
ISM MIR polarisation p rofile presented in lWright et alj l2002h and 
IWright & Glassel l2005h . This provided very similar fits to the data. 

In summary, both the conventional and polarisation spectra re¬ 
quire the presence of a material with optical properties much like 
Sic, inclusive of a high band strength, a central wavelength near 
11.3 pm and full width of a few microns. We thus assess the car¬ 
rier of the additional absorption from approximately 10 to 12 pm 
in SVS13 to be SiC itself. But given the variation in the optical 
properties of specific types of SiC found in the literature a reliable 
estimate of its abundance is difficult. This is apart from a general 
statement that it is higher than in any other young star envelope or 
disc of which we know. 


4.1.3 Crystalline silicates 

The ISO crystalline silicate revolution ilager et alj Il998l) came 
about because the ISO satellite was equipped with instruments ca¬ 
pable of making measurements at longer wavelengths, where crys¬ 
talline silicate features are not overwhelmed by those of more abun¬ 
dant and w armer amorphous s ilicates, as is usually the case in 
the Wband l lMolster et alj|201^ . There are three plausible mecha¬ 
nisms adept at producing crystalline silicates. 


1. Evaporation and re-condensation Under high temperature 
and pressure condition in the immediate vicinity of the central star 
(< 1 au), du st grains can be evaporated and then re-condense as 
crystals (e.g. jGrossma 

2. Thermal annealing There have been at least three probable 
heat sources suggested in cir cumstellar d iscs [accretion luminosity 
in the inner region ( ~1 au ), iGaill l l2004h : shock heating at a few 
au, iHarker & Deschl | |2002|): and disc surfac e layer annealing dur¬ 
ing accretion outburst, Abraham et al.l ( l2009l) 1. Whatever the actual 
heat source(s) may be, the essence of the mechanism is simple an¬ 
nealing (i.e., heating of amorphous silicates and subsequent cooling 
to re-order the lattice structure). 

3. Low-temperatu re crystallisation Put forwa rd by Yamamoto 
and co-workers (e.g.. lYamamoto & ChigafcoOSh . provided that an 
amorphous silic ate core is coated with an organic (carbonaceous) 
refractory layer l lKimura. Mann & Jessbergen2003h . a chemical re¬ 
action can be triggered at a low temperature (a few hundred K) to 
re-arrange the silicate lattice structure. 

It is interesting to note that, so far, most detections of crys¬ 
talline sili cate features towards young stars have been in emis- 
sion (e.g., IWatson ^ l2009l : lOlofsson et al.l l2009l : ljuhasz et al.l 


I 2 OIC ), with only one recent r^ort of clear detection of absorp¬ 
tion I Poteet et al. II 20 111) to datejj a fact that may or may not favour 
one crystallisation mechanism over the others. Inclusions of the two 
most abundant crystalline silicate species, forsterite (Mg 2 Si 04 ) and 
enstatite (MgSiOs), in absorption, substantially improve the model 
fitting in the current study. Equilibrium calculations suitable for the 
inner disc regions (< 1 au. lGai I 2 OO 4 show that these magnesium- 
rich crystalline silicates should dominate. 

However, the crystalline silicate features are observed in ab¬ 
sorption here and therefore the dust grains responsible for the sig¬ 
natures are likely to reside in the cold, outer regions. Apart from 
the low-temperature crystallisation described above, some means 
of heating is required to turn the amorphous lattice structure into 
crystalline, and if heating is only possible in the vicinity of the 
central star, an effective transport/mixing mechanism is required 
for crystalline silicates to be found in the remote, cooler regions. 
Indeed, crystalline silicates appear to be concentrated in the cen¬ 
tre of many T-Tau syst ems, with more tha n 90 per cent of the 
84 objects observed by I Watson et alj l l2009tl showing strong 10- 
qm crystalline silicate features but only about 50 per cent hav¬ 
ing detectable 20-38-qm signatures, most probably indicating the 
temper ature range of the c r ystalline silicates (T ~ 200 K at < a 
few au. lWatson et alJl2009l) . foiofsson et ^ l l2009l) . who, as part of 
the Spitzer cores-to-discs legacy programme, observed 108 young 
stars (66 of which were known T-Tau stars), identified what they 
termed the ‘crystallinity paradox’, wherein they found the cold 
crystalline features at wavelengths > 20 qm more frequently than 
those at ~10 qm. However, when they re-examined it by analysing 
a further 58 Spitzer IRS spectra, they found a simultaneous en- 
hancem ent of crystallinity in both the inner, warm a nd outer, cold 
regions l lOlofsson et alJl201^ . ISargent et alj i2009bh also reported 
comparable crystalline silicate abundances in the inner and outer 
disc regions in a Spitzer IRS study of 65 T-Tau stars. Although 
details may differ slightly, an undisputed universal outcome of all 
these investigations is the existence of processed (i.e., crystalline) 
materials in the outer, cooler parts of the circumstellar environment, 
as well as in the warmer, inner regions in the vicinity of the central 
star. 

Comets contain crystalline silicates (e.g.. Comet Hale-Bopp, 
ICrovisier et al ][l99i; 1997). The most definitive confirmation was 
once again found in the Stardust sample return from Comet 
81 P/Wild 2, the crystall ine silicate mass frac tion of which was 
as high as ~0.5-0.65 l lWestphal et alj l2009l) . The 81 P/Wild 2 
speci mens also included chondrule-like particles jNakamura et al.l 
l2008h . a clear indication that materials that constitute the 
comet experienced high temperatures in the past. Another clue 
to the apparent heating episode(s) in the Stardust samples is 
the pre sence of Calcium-, A luminium-rich Inclusion (CAI)-like 
grains jZolenskv et al.l 1^0061) . However, as mentioned in 94.1.21 
81P/Wild 2 belongs to the short-period JFCs, which are believed 
to originate in the Kuiper Belt, extending from the orbit of Nep¬ 
tune (at ~30 au) to approximately 50 au, in which the surface tem- 
perature of the most famous of i ts resident, Pluto, is only ~40 K 
(IStern. Weintraub & Festoulll993h . 


^ A tentative identification has been made for d eeply embedded low -mass 
protostars IRAS 03445+3242 and H H46 IR S iBoogert et al.| l2004ll . and 
for the high-mass YSO AFGL2591 iAitken et al.lll988ll . Crystalline sili¬ 
cate absoiption fe atures towards oth er types of objects, e.g., ultra-luminous 
infrared galaxies iSpoon et al.ll2006h . have also been repoifed. 
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To overcome these facts at odds with each other (existence 
of high-temperature phases in cold regions), a number of mecha¬ 
nisms have been proposed to mix/transport radially the dust grains 
transformed near the cent ral star. These include: t he bipolar out¬ 
flow ‘X-wind’ m odel (e.g., ShUj_Shan^_& Led 19961). accret ion disc 
turbulence (e.g., iGaillI 2 OOII : iBockelee-Morvan et al.l l2002h. spira l 
arms in marginally gravitationally unstable discs (e.^ Bosj2004 ), 
and outward transp ort in the disc mid-plane (e.g.. iKeller & Gai] 
|2004 lOeslal 120071) . Now, both chondrules and CAIs are sizable 
entities that often measure hundreds of pm and so me as large as 
mm (chondrules) or even cm (C AIs) in diameter (e.g.. lScott & Kroli 
l2005l : lMacPherson et al.l20^. Although particles that large in th e 

81 P/Wild 2 capture are rare fflorz et al.l2006l : lBurchell et al.l2008h . 

analogous to chondritic meteorites that contain large chondrules 
and CAIs, a similar size distribution may also be expected to ex¬ 
ist in their probable nat al bodies such as comets l Meieill2oI4) . A 
recent modelling effort ( IHughes & Armitagell201C ) has found that 
particles of the size of chondrules and CAIs (a few mm or larger) 
cannot easily be transported out to the comet-forming region, re¬ 
gardless of their initial location. If this is the case, an in-situ pro¬ 
duction of melted grain may have to be invoked. We will come 
back to this point in i l4.2l 


4.1.4 Annealed Si02 


Aforementioned equilibrium calculations also show that, when 
the Mg/Si abundance ratio is < 1 (or perhaps after a signifi¬ 
cant amount of Mg has been taken up by fo rming forsterite), 
Si02 ( and enstatite) becom es a major condensate jperrarotti & Gaill 
l200ll) . lFabian et al] ( l200(]|) found that when they annealed enstatite 
(MgSiOs) smoke at 1,000 K for 30 hours, the end product consisted 
of crystalline forsterite (Mg 2 Si 04 ), tridymite (crystalline Si 02 ) and 
amorphous silica (S i O?). T his is a case similar to that reported by 
ISowen & Anders^ ( Il914t) . who cooled a mixture of composition 
MgSiOs from the liquid state and found that, first forsterite sepa¬ 
rated out and, after further cooling, a mixture of silica and enstatite 
condensed. As crystalline silicates also exist, it may be deduced 
that some form of heating has melted MgSi 03 and decomposed it 
into Mg 2 Si 04 and Si02. We further speculate that, once again, tak¬ 
ing the presence of crystalline forsterite and enstatite as evidence, 
the cooling rate is just so that it would allow re-ordering of lat¬ 
tice structures to produce annealed Si02. Or, that crystalline Si02 
is dominant over its amorphous counterpart in the stages immedi¬ 
ately after their formation. 

As mentioned earlier, SiO? has been ide ntified in the MIR 
spectra of some young stars. iHonda et alj ( l2003h observed a T-Tauri 
star Hen 3-600A in the A-band and their be st-fit model included 
tt-quartz (i.e., a type of cry stalline Si02, from lSoitzer & Kleinmanl 
1961^ Sargentetaljj 2006h us ed optical constants of ff-quartz from 


Wenrich & ChristensenI ( Il996l) to analyse A-band spectra of 12 T- 


Tau stars. They found a small amount (up to a few per cent by mass) 
of crystalline Si02 in about a half of the objects studied. Amor- 
phous silica i Henning & Mutschkel Il997h and sil ica-rich glass 
I Koite et ^ 1 19891) emissivities were ruled out by I Sargent et al.l 
l 2009ah . his time observi ng 5 more, differen t T-Tau stars, but in¬ 
clusion of annealed Si02 dpabian et al.ll2000h gave a better fit. A 
Spitzer IRS spectrum of a /1-Pictoris analogous (i.e., debris-disc) 


® Note that grains in question here are not crystalline silicates so the low- 
temperature crystallisation of Yamamoto and co-workers mentioned earlier 
would not apply. 


star HD172555, which suggested t he presence of am orphous sil¬ 
ica and SiO gas, was examined by iLisse et al.l ( l2009h . The 12.6- 
pm feature that would have indicated the existence of crystalline 
silica was not detected. They proposed planetesimal-scale (~km- 
size) hyper-velocity (> 10 km s“*) impacts as the mechanism 
that pr oduced fine (amorpho us) silica dust and SiO gas. More re¬ 
cently, iFuiiwara et al.l l l2012h obtained a Spitzer IRS spectrum of 
another debris-disc star HD 15407A and concluded that incorpo- 
rating an almost e qual amount of fused (i.e., amorphous) q uartz 
dKoike et alJll989h and annealed silica dPabian et aljl2000h pro¬ 
vided the most satisfac tory fit. They also favoured a similar sce¬ 
nario to that adopted bv iLisse et alj d2009t) for HD 172555. 

Although the number of samples is rather limited and most 
studies mentioned above did not set out specifically to probe crys¬ 
talline or otherwise of silica dust structure, overall it is tempting 
to conclude that younger stars (T-Tauri stars, < 10 Myr) are more 
likely to possess more crystalline than amorphous Si 02 , while with 
slightly more evolved sources (debris-disc bearing stars, > 10 Myr), 
the balance is tipped in the other direction. Also, if hyper-velocity 
planetesimal impacts are indeed required for the production of 
amorphous silica, this alone points to the latter stages in the evolu¬ 
tion of protoplanetary discs. If this scenario is correct, then the nat¬ 
ural implication may be that the circumstellar environment of (still 
deeply embedded and therefore much younger) PMS star SVS13 
should contain mostly crystalline Si 02 but not much (if at all) of its 
amorphous sibling. 


4.2 Possible crystallisation mechanism and planet formation 

We have demonstrated that crystalline silicates (forsterite and en¬ 
statite) and annealed Si 02 exist in the cold, outer regions of cir¬ 
cumstellar environment of the PMS star SVS13. Although we did 
not directly detect larger chondrules and CAIs in the current study, 
they may exist if the typical dust size distribution in chondritic 
meteorites also holds in the remote, cooler regions (at least in the 
comet-forming region at a few tens of au from the central star). If 
that is the case, an in-situ transformation mechanism of dust grains 
is probably required. 

The runaway growth of planete simals is the ‘standard ’ model 
for terrestrial planet formation (e.g.. [Raymond et^l2014l) . Outer 
gas giants, like Jupiter and Saturn, can be formed in a similar 
fashion via what is called the core-accretion (CA) model (e.g., 
iHelled et al.ll2014l) . But also feasible is if the remote, cold regions 
of the protoplanetary disc become gravitationally unstable and frag¬ 
ment, forming self-gravitating clum ps, some of which will ulti¬ 
mately become gas giants (see e.g., iHelled et "^iTOldl and refer¬ 
ences therein). 

According to the latter mechanism, commonly referred 
to as the disc instability (DI) model, fragmentation can take 
place at > 50 au and clumps ranging from a few to ~10 Mj 
(Jupiter masses) may form (e.g. , iBolev et^l2010l : iForgan & Ricel 
I 2 OIIL [Rogers & WadslevI [2012[) . As these massive clumps con¬ 
tract and migrate inwards, the internal temperature rises and it 
can reach > 1,000 K, sufficient to produce crystalline silicates, 
chondrules and CAIs (and dissolving amorphous carbon coat¬ 
ing on Sic grains) in situ. When they are tidally disrupted, 
these thermally altered particles may be distri buted at a wide 
range of radii to be re-captured by, e.g., comets |[Vorobvov[[201 ll : 
[Navakshin. Cha & BridgedPTOllI) . 

Of course, thermal processing of dust species can still occur 
in the immediate vicinity of the central star and the dust grains 
may be transported out to the comet-forming region of the pro- 
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toplanetary disc. One of the plausible dust heating processes that 
could produce chondrules and CAIs (and at the same time a pos¬ 
sible transport/mixing mechanism), is the X-wind model of Shu 
et al. (i.e., jets near th e central star). But this has been criticised 
by IPesch et"^ ( l2010h. who prefer shock-h eating, d ue to PI but 
at a f ew au (e.e.. IPesch & Connollvl [ 2 OO 2 I : see also IPesch et al.l 
l2012h . However, as mentioned in 34.1.31 it is difficult, if not impos¬ 
sible, to move large particles such as chondrules and CAIs out to 
the c omet-forming region, where ver their initial location may have 
been iHughes & Armita^l2010h . Even if shocks can be triggered, 
the gas density at the distance of the Kuiper Belt (> 30 au) is such 
that chondrule formation by s hock-heating might not be feasible. 
Ilida. Nakamoto & Susal 1200 ih estimate the pre-shock density of at 
least lO*'^'^ cm“^ is required , while both the minimum-mass So¬ 
lar Nebula mod el of iHayashil ll9Mb and the minimum-mass extra- 
Solar nebula of IChiang & LaughlinI l l2013h predict a much lower 
density of the order ~10*®“*' cm“^ at 30 au. 

We recall that SVS13 is a close binary, with a separation of 
~0.3 arcsec » 71 au. Theory and model simulations depict some 
common features in the early stages of binary evolution: three sep¬ 
arate discs (individual circumstellar discs around each component 
of the binary system and a circumbi nary disc) and a clearing be¬ 
tween the two classes of discs (e.g., Art ymowicz & Lubowlll994l : 
iBate & Bonnel]|[l997l : [Giinther & Klevll2002l) . Under certain con¬ 
ditions, however, one binary component (primary) will gain a large 
circumstellar disc, while the other (secondary) onl y a small one (or 
none at all), with little to no circumbinary disc jBate & Bonnelll 
Il997h . The primary is usually the do minant dust continu um source 
in PMS binary/multiple systems (e.g.. lHarris et alj2012h . and apart 
from a few outstanding cases, circumbinar y discs appear to be elu¬ 
sive jMonin et alJl200'4 : lHarris et al.ll2012t) . In the case of SVSI3, 
only one mem ber of the binary (V LA 4B) exhibits significant 
dust emission jAnglada et alj 120041) and just a single object has 
been found in the NIR (which Hodapp & Chinill2014l identify as 
VLA 4B) as we do here in the MIR. Furthermore, the unresolved 
radio source exhibits a larger flux density than the sum of the re¬ 
solved individual bi nary components, im plying the existence of an 
extended envelope jAnglada et alj[200^ . Also commonly seen in 
snapshots from simulations are spiral arms in all three discs due to 
the tidal interactions of the central binary. Whether such gravita¬ 
tional interplay and structures i nhibit or enhance Pl-induced frag¬ 
mentation is a matter of debate jMaver. Boss & Nelsonll20l'^ . 

The ‘clearing’ or a gap between the circumbinary and the cen¬ 
tral circumstellar discs is also a product of the tidal interactions 
between the binary components, as the protoplanetary discs are 
stripped of materials and truncated to sub stantially smaller sizes 
comp ared to th ose around single star s (e.g.. lArtvmowicz & Lubowl 
Il994h . Indeed, lAnglada et al.l j2004l) found a compact structure at 
the position of VLA 4B of radius ~30 au, clearly smaller than 
that often measured for discs around single stars of > 100 au (see 
e.g.. I Williams & Ciezall20ll] and references therein). Discs are not 
only truncated but, in such a dynamic environment, they are also 
rapidly dispersed within a short period of time (< 1 Myr), con¬ 
siderably shorter than a characteristic di sc lifetime of > 2.5 Myr 
normally expected for single stars (e.g.. lMamaiek|[200^ . though 
the effect could be slightly milder for a binary separati on of 
~70 au for SVS13 jMonin et alj l2007l : ICieza et al.l 1200^ . Al¬ 
though this type of disruptive impact might be thought to affect the 
planet formation negatively, Bonavi ta & Desidera (2007; see also 
iBonavita. Desidera & GrattorJ 120101) concluded that the detection 
frequencies of planets in sing le and bin a ry sys tems are statistically 
identical. Most interestingly, IPuchenel j201(]h showed that, while 


binaries with separation wider than 100 au are indistinguishable 
from single stars in terms of trends regarding their protoplanetary 
discs, debris discs, and planets, close binaries in the range from 
5 to 100 au exclusively host giant planets (i.e., those with mass 
> 1 Mj), or exhibit a distinct lack of planets smaller than 1 Mj. 
The two modes of planet formation discussed earlier present no- 
tably different time scales: CA with > a few Myr and DI < 1 Myr 
jHelled etalj2014h . Considering the short lifetime of circumstellar 
discs in close binaries, it could be the case that DI is the only ef- 
fective means o f forming planets in them, and if so, the finding of 
IPuchenel j2010h might just be a natural consequence of this. 

Recen t observations o f the young low-mass binary system 
GG Tau A jBeck et al.l2012t) have revealed that a ‘streamer’ of dust 
and gas, originating from the circumbinary disc and feeding the 
inner circumstellar discs can shock-excite H 2 gas at the interface, 
where the gas temperature was found to be in excess of 1,500 K. 
Unlike the shock-heating of particles (due to DI) rejected earlier 
because of the low density in the outer parts of the circumstellar 
disc, we surmise that a stream of material from the circumbinary 
disc piling up at the outer rim of the protoplanetary disc, combined 
with a higher shock velocity generated b y streamers (20-30 km s~’ 
for streamers vs. < 1 ( ) km s~‘ for DI, IPesch & Connollvl 1 20021 : 

iHarker & Desch|[2003 : lBeck et alj|2012h . could make this a viable 

mechanism for melting/annealing dust grains in situ in the remote 
regions of circumstellar discs. 

We therefore suggest that, amongst the processes that are ca¬ 
pable of thermally altering dust grains, an in-situ mechanism is 
favoured for SVS13. Furthermore, the DI model for giant planet 
formation - provided that fragmentation occurs in circumstellar 
discs in bin ary systems - alon g with the shock-heating mechanism 
detected bv iBeck et alj j2012h . is the leading candidate for produc¬ 
ing crystalline silicates (and chondrules & CAIs), annealing Si02, 
as well as possibly dissolving the feature-concealing amorphous 
carbon layer from the SiC core. It remains to be seen if the scenar¬ 
ios just presented are rare or only occur in a short-lived evolutionary 
phase in shaping the distinctive A-band spectrum of SVS13. To this 
end, we are embarking on the search for similar objects and at the 
same time, it is hoped that work both in the laboratories and theo¬ 
retical studies will advance sufficiently in the near future to solve, 
for example, the missing SiC problem. 


5 CONCLUSIONS 

The A-band spectrum of the low-mass PMS close binary system 
SVS13 is presented. The unique and complex spectrum is best 
modelled by a mixture of amorphous silicates, crystalline forsterite, 
crystalline enstatite, annealed SiOi, and, most intriguingly of all. 
Sic, which has never before been identified in the circumstellar en¬ 
vironment of a young star. The requirement for the inclusion of sil¬ 
icon carbide is also affirmed from the MIR polarimetry perspective. 
All these signatures are seen in absorption, implying the dust parti¬ 
cles reside in remote, cold regions of the circumstellar environment. 
Speculation is made on possible origins of the dust species, espe¬ 
cially those that have been processed and altered thermally, namely 
crystalline silicates and annealed Si02. The Dl-induced fragmen¬ 
tation and the subsequent contraction and disruption of the giant 
planet embryos, along with the newly discovered shock-heating 
mechanism at the interface between the circumbinary and circum¬ 
stellar discs, may be able to modify those dust grains in situ in the 
outer, cooler parts of the circumstellar environment. The heating 
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episode(s) provided by one or both of these processes can also fea¬ 
sibly reveal the hidden SiC feature. 
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